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Abstract 

In this paper we propose a transform method to compute the prices and greeks 
of barrier options driven by a class of Levy processes. We derive analytical ex- 
pressions for the Laplace transforms in time of the prices and sensitivities of single 
barrier options in an exponential Levy model with hyper-exponential jumps. In- 
version of these single Laplace transform yields rapid, accurate results. These 
results are employed to construct an approximation of the prices and sensitivities 
of barrier options in exponential generalised hyper- exponential (GHE) Levy mod- 
els. The latter class includes many of the Levy models employed in quantitative 
finance such as the variance gamma (VG), KoBoL, generalised hyperbolic, and 
the normal inverse Gaussian (NIG) models. Convergence of the approximating 
prices and sensitivities is proved. To provide a numerical illustration, this trans- 
form approach is compared with Monte Carlo simulation in the cases that the 
driving process is a VG and a NIG Levy process. Parameters are calibrated to 
Stoxx50E call options. 
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1 Introduction 



Barrier options are derivatives with a pay-off that depends on whether a reference 
entity has crossed a certain boundary. Common examples are the knock-in and knock- 
out call and put options that are activated or des-activated when the underlying crosses 
a specified barrier-level. Barrier and barrier-type options belong to the most widely 
traded exotic options in the financial markets. Whereas knowledge of the marginal 
risk-neutral distribution of the underlying at maturity suffices to obtain arbitrage 
free prices of European call and put options, the valuation of barrier options requires 
specification of the risk-neutral law of the underlying price process, as it depends on 
the first-passage distribution of this process. At least as important as the calculation 
of prices is the evaluation of the sensitivities of the prices to various model-parameters 
(the greeks) for which the law of the process is also required. 

A class of models that has been shown to be capable of generating a good fit of 
observed call and put option price data is formed by the infinite activity Levy mod- 
els, such as extended Koponen or KoBoL pjj, variance gamma [29J, normal inverse 
Gaussian [6] and generalised hyperbolic processes [H] . This class of models has been 
extensively studied and we refer for background and further references to the books by 
Boyarchenko and Levendorskii [T3], Cont and Tankov p!B| and Schoutens In this 
paper we consider barrier options driven by Levy processes with a completely mono- 
tone Levy densities on each half-axis (which we call generalised hyper- exponential) . 
This class contains many of the Levy models used in financial modelling as the fore- 
mentioned ones and also jump-diffusions with double-exponential jumps, as shown in 
Section [2l 

The calculation of first-passage distributions and barrier option prices in (specific) 
Levy models has been investigated in a number of papers. Geman and Yor [22] calcu- 
lated prices and deltas of double barrier options under the Black-Scholes model. For 
spectrally one-sided Levy processes with a Gaussian component Rogers [Slj derived a 
method to evaluate first-passage distributions. Kou and Wang [26], Lipton [28] and 
Sepp [35| followed a transform approach to obtain barrier prices for a jump-diffusion 
with exponential jumps. In the setting of infinite activity Levy processes with jumps in 
two directions Cont and Voltchkova p21 investigated discretisation of the associated 
integro-differential equations; Boyarchenko and Levendorskii [12] employed Fourier 
methods to investigate barrier option prices for Levy processes of regular exponen- 
tial type; Asmussen et al. |5j priced an equity default swap under a CGMY model, 
by fitting a hyper-exponential density to the CGMY Levy density. In this paper we 
adopt the latter approach. As first step we obtain, in a Levy model with hyper- 
exponential jumps, analytical formulas for the Laplace transform in time of knock-in 
and knock-out barrier option prices by exploiting the availability of an explicit Wiener- 
Hopf factorisation. Using these results we also establish analytical formulas for the 
corresponding sensitivities with respect to the initial price (delta and gamma) and 
the time of maturity (theta) up to one Laplace transform in time. The actual bar- 
rier option prices and sensitivities are then obtained by numerically inverting these 
single Laplace transforms, using Abate and White's algorithm [1], yielding fast and 
accurate results. Since hyper-exponential Levy processes are dense in the class of 
generalised hyper-exponential Levy processes, the idea is to approximate the barrier 
option prices and sensitivities in a generalised hyper-exponential model by those in 
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an appropriately chosen hyper-exponential Levy model. At this point it is worth re- 
marking that for a general Levy process the Wiener-Hopf factors are not available in 
analytically tractable form (as they are expressed in terms of the one-dimensional dis- 
tributions that are generally not available) and, furthermore, even if the Wiener-Hopf 
factors have been obtained still a three-dimensional Laplace/Fourier inversion would 
be needed to obtain the knock-in and out call option prices (see e.g. jlBl p.372]). 

Following the approach described in the previous paragraph barrier option prices 
and sensitivities under a generalised hyper-exponential Levy model can, at least in 
principle, be approximated arbitrarily closely. Indeed, we will prove that, when a 
sequence of hyper-exponential Levy processes weakly converge to a given generalised 
hyper-exponential Levy process of infinite activity, the corresponding first-passage 
times converge in distribution, and the barrier option prices and smoothed sensitivities 
converge pointwise to those of the limiting model. We illustrated this approximation 
procedure by implementing it for the exponential variance gamma (VG) and normal in- 
verse Gaussian (NIG) models, with parameters calibrated to StoxxSOE options. Using 
a least-squares algorithm to minimize the root mean square error of the approximating 
density with respect to the target density over a grid, using 7 upward and downward 
phases, we determined the parameters of the approximating hyper-exponential Levy 
densities; the resulting hyper-exponential Levy processes we employed as approxima- 
tions to the VG and NIG processes. We calculated the prices and sensitivities of 
European and American digital options and down-and-out put options following this 
approach, and also by Monte Carlo simulation. We found a general agreement be- 
tween the results, with relative errors in the range of 0.5-2.5% for prices and deltas, 
some distance away from the barrier. Numerical experiments showed that closer to the 
barrier errors may be larger, especially for the delta, suggesting that a more accurate 
approximation of the target density by a hyperexponential density would be needed 
to reduce the size of the error, which could be achieved by increasing the number of 
phases or by employing one of the methods from the approximation theory of real val- 
ued functions. The phenomenon of larger errors in the vicinity of the barrier was also 
observed by Kudryatsev and Kudryatsev and Levendorskii [2^ when approximating 
first touch digital option prices under a NIG Levy model using the Kou model. It 
would be desirable to analyze the dependence of the error on the different parameters 
and the distance to the barrier, and how the presented approach compares to alterna- 
tive approaches such as finite difference schemes, but, in the interest of brevity, these 
questions are left for future research. 

The remainder of the paper is organized as follows. In Section [2] the 'generalised 
hyper-exponential Levy model' is defined and it is shown that many of the existing 
Levy models used in quantitative finance are contained in this class. Results on the 
Wiener-Hopf factorisation and first passage times for processes from this class can be 
found in Section [31 Analytical results for the Laplace transforms in time of barrier 
option prices are obtained in Section HI Using these results explicit expressions are 
derived in Section for the Laplace transforms of the theta, delta and gamma of the 
barrier an digital options. In Section [6] numerical results are presented for prices and 
sensitivities of barrier options in a VG and a NIG model respectively, with convergence 
results presented in Section [71 Some proofs are relegated to the Appendix. 
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2 Model 



Let {X{t),t > 0} be a Levy process, that is, a stochastic process with independent 
and homogeneous increments, normahsed such that X{0) = 0, defined on some ap- 
propriate probabihty space {Vt,J-,P). For background on the fluctuation theory of 
Levy processes and the apphcation of Levy processes in quantitative finance we re- 
fer to Bingham [9] and Cont and Tankov [T6], respectively — Sato [M] and Bertoin 
[H] are general treatments of the theory of Levy processes. Assume that X satisfies 
E[%'^^^'^] = e^''"'^^* for all t > 0, where r and d are non-negative constants representing 
the risk-free rate of return and the dividend yield, and consider the model for the 
risk-neutral evolution of the stock price S given by 

= 5oe^W. 

As a consequence of the independent increments property of X, e~^^~''^^^ S{t) is a 
martingale (under the measure P and with respect to its natural filtration). We will 
restrict ourselves to the following class of Levy processes: 

Definition. A Levy process is said to be generalised hyper- exponential (and we will 
denote this class of processs by GHE) if its Levy measure admits a density k of 
the form k{x) = A;_|_(a;)l|2,>o} + fe-(^a^)l{x<o} where k+, A;_ are completely monotone 
functions on (0, co) and 1a denotes the indicator of a set A. 

In view of Bernstein's theorem a Levy process X is a member of the class GHE if and 
only if its Levy density k is of the form 

fe(x) = l{,>o}/ e-">+(dn) + l|,<o} / e-M^_(d^.) (1) 

for some measures fi+,fi- on (0, oo) and (— oo,0) respectively. The mass of the inter- 
val [a,b], under the measure /i+ corresponds to the frequency of positive exponential 
jumps of mean sizes between 1/a and 1/6. A similar statement holds true for the 
negative jumps and In the case that the measure /i+ is a convex combination of 
point-masses, the locations Ui > and sizes pi > of the point-masses respectively 
correspond to average size and frequency pi of the exponential jumps. Below we 
show that many of the Levy processes employed in financial modelling are generalised 
hyper-exponential by explicitly calculating the corresponding measures ii±. In partic- 
ular, the class GHE contains the time changes of Brownian motion by a generalised 
hyper-exponential subordinator. 

Proposition 2.1 Let W be a Brownian motion and Y an independent generalised 
hyper- exponential subordinator. Then, for /i G M, X{t) = W{Y{t)) + f^Y{t) is a Levy 
process in the class GHE. 

Proof Denoting by p the Levy density of Y and /Uy the measure in the representation 
([1]), Theorem 30.1 in Sato (1999) implies that the Levy density k of X is given by 

k{x) = / _^e-(^-'^^)'/(2.)p(^)ds. 
Jo v27rs 



3 



An interchange of the order of integration yields that 



OO 1 fOO 



k{x) = I — ^e-(^-^^) / e-""/iy(du)ds 
/o v27rs Jo 

f OO /"OO 



/"OO /"OO 1 

/ / ^=e-(^-'^'^)'/(2.)g-s«^^^^(^^) 
Jo Jo v27rs 

OO 1 ^ 



/O ^J jj? + 2U 

and the claim follows. □ 

The following result gives necessary and sufficient conditions on the measures 
such that A; in ([T]) is a Levy density: 

Proposition 2.2 Let ^ he a measure on M\{0} and set //±(dx) = l{-|-j.>o}A*(dx). 
Then k defined in ^ is a Levy density if and only if 

-i- A fiidu) < OO. (2) 

Proof By interchanging the order of integration it can be verified that 

OO /"OO /"OO Q — u 

e~'^^ fi{du)dx = / M(dn), 

./n u 



1 /"OO 

2 / „—ux 



x' I e-"Xdn)dx = / ^- — (l + l + \]^{du). 



Further we note that Jq^ e "^^(dti)dx is bounded below and above by ^^^^ and 
^^^^^ respectively. In view of these calculations we read off that A: in ([1]) satisfies the 
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integrability condition /[I A x'^]k{x)dx < oo (which is the requirement for to be a 
Levy density) if and only if ([2]) holds. □ 

In the following examples we explicitly determine the measures 

Examples. We shall write A;+ for the density k restricted to the positive half- axis. 

• Double exponential model (Kou [25]) For a jump-diffusion model where the 
jumps follow a double exponential distribution, /i+ is a point-mass located at the 
reciprocal of the mean jump sizes: 

k+{x) = A+a+e"""^^, Ai+(dtt) = 5a+{du), 

where , > 0. 

• Hyper-exponential model (e.g. [3]) This model is an extension of Kou's model 
by allowing for exponential jumps with a finite number of different means. For pf , af, 
and > with Yll^=i = 1 we thus get 

n n 

k+{x) = A+ 5]p+a+e-"^^ t,+ {du) = A+ ^p+a+5„+(dn). 

j=l i=l 
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• KoBoL/CGMY model ( [IH [H] ) In the KoBoL model (also called CGMY model), 
the measure fi has a continuous density k, and it holds that 

where C, M, y > 0. The form of invokes the definition of the gamma function 

1 r°° 



x^+i Jo r(l + Y) 

In particular, for the variance gamma model, we set Y = and get k-^-{x) = Cx~^e~^^ 
and fi+{du) = Cl{„>Af}dtt. 

• Meixner model (e.g. [36]) The Levy density of a Meixner Levy process is given 
by 

xsmhivrx/a) \x\ 

^ ' ' n=0 ' ' 

where 5, a > 0, — tt < /3 < tt, so that we find that 

oo 

^+(du) =25^ l{«>{27rn+7r-/3)/a}du. 
n=0 

• Normal-inverse Gaussian (NIG) ( Bar ndorff Nielsen [6]) In the NIG model, the 
measure [i has a density which reads as 



k+ix) = e^" /i+ dn = + /? /a 2 - 1 i/2i ^u, 

TT X IT I - /-J 



where a>/?>0, 5, C>0 and Ki is the McDonald function. The form of is based 
on the following representation (see [2]) of Ki 

Ki{x)=x / e-^'''(i;2 - 1)1/2^^;. 



• Generalised hyperbolic (GH) (Eberlein et al. [H]) The GH process can be 
respresented as time change of Brownian motion by a generalised inverse Gaussian 
(GIG) subordinator. The Levy density of GIG is a generalised gamma convolution 
which means in particular that it is of the form ([T]). 



3 First passage theory 

First passage distributions are an essential element in the valuation of barrier options. 
In this section we briefly review elements of the first passage theory for Levy processes 
that will be needed in the sequel. 
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3.1 Wiener-Hopf factorisation 

The distributions of the running supremum and the running infimum of X are hnked to 
the distribution of X{t) via the famed Wiener-Hopf factorisation of the Levy exponent 
of X, denote by ^{u) = logS[e'"^(^)]. For v e C, let Q{v) denote the imaginary part 
of u. 

Definition. A Wiener-Hopf factorisation of ^ is a pair of functions F^,F^ that 
satisfies, for u E M and q > 0, the relation 

q{q-^{u))-'=F+{u)F-{u), (3) 

where, for fixed q > 0, n i— > F^{u) are continuous and non-vanishing on ib$5(u) > 
and analytic on '^{u) > with F^{0) = 1, and 



q 

lim 1 

<? 



lim F^(u) = 1. (4) 



Denote by X[t) and Xif) the running supremum and infimum of X up to time t, 
X{t) = supX(s), Xit) = inf 

s<t 

and let 5"^$^, q~^(^~ be the respective Fourier-Laplace transforms 

$+(n) = / ge-«*^[e™^W]dt, ^^{u) = / ge-«*^[e'"^W]dt. 

It is not hard to verify that satisfy the condition Rogozin [33J established the 
following result: 

Theorem 1 (Rogozin) <I>^ (u) , (ti) is the unique Wiener-Hopf factorisation of^. 

If the Levy measure has support in both half-axis the factorisation is generally not 
explicitly known, except if the Levy measure of X is of a particular form. The obser- 
vation that the Wiener-Hopf factorisation and related first passage distributions are 
tractable for a random walk in the case that the jump distribution on the positive 
half-axis follows an exponential distribution can already be found in Feller [20]. See 
also the books Borovkov [10] and Asmussen f^ for background on the Wiener-Hopf 
factorisation for random walks. Mordecki |30j and Lipton j28j calculated the Laplace 
transforms of first passage time distributions for a Levy process with hyper-exponential 
jumps. Asmussen et al. [3] derived explicitly the Wiener-Hopf factorisation if the jump 
distributions of the Levy process are of phase-type. 

In the sequel we will employ the known factorisation results for a jump-diffusion 
with Levy density (also to be called a hyper- exponential jump- diffusion (HEJDJ) given 
by 



k{x) = X+Y^pfafe-'^^ n{,>o} + A- ^pTaTe-«"-l|,<o}, (5) 
i=i j=i 

where af, X^pf > with Y17=iPf ~ "^^^ corresponding Levy exponent then reads 



as 



*(n) = /.ni - ^u' + A+ Y^pt -l)+X- ^Pj i " 1 I • (6) 

2 ^ \al — u\ J ^ ■' \a.- j 
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The function \I' in ((61) can be analytically extended to the complex plane except the 
finite set {— iQ^,ia~,i = 1, . . . ,n^}, and this extension will also be denoted by 
Denoting by pf = pf{q), i = 1, . . . , m+ and p~ = p~{q), j = 1, . . . , m~ the roots of 
'!'(— is) = q with positive and negative real parts respectively, the Wiener-Hopf factors 
and $~ are explicitly given as follows: 



Ul 



'^3 



Note that as a consequence of the intermediate value theorem and the fact that s i— > 
\1'(— is) is a rational function with real coefficients and single poles, it follows that the 
roots of the equation ^'(— is) = q are all real and distinct if g G M+. Performing a 
partial fraction decomposition and termwise inverting the terms yields that the time- 
Laplace transforms of the distributions of X{t) and X(t) are given b}0 



J e-^'P(X{t) < z)dt = - I 1 - A 



poo 1 / 

e-'^*P(-X(t) < z)dt =-U-J^A-^ 
where the coefficients Af and A~ are given by 



2 > 0, 



z > 0, 



At = ^ ^^t-^, A-, 



n:::i(i+^ 



ni)=i,D5^i ( 1 



(7) 



3.2 Generalised hyper- exponential Levy processes 

For a generalised hyper-exponential Levy process X the Wiener-Hopf factorisation is in 
general not known explicitly. However, due to its special form ([T]), its Levy density can 
be approximated arbitrarily closely by a hyper-exponential density (by approximating 
the measures p± by sums of point-masses), and for the resulting hyper-exponential 
Levy process the Wiener-Hopf factorisation is explicitly given in the previous section. 
As hyper-exponential Levy processes are 'dense' in the class GHE, the Wiener-Hopf 
factorisation of X can be approximated arbitrarily accurately using this idea. 

More formally, a sequence of approximating processes (X("'))„ can be explicitly 
constructed as follows. Denote by k and the Levy density and the diffusion coeffi- 
cient of X, respectively. Fix a sequence {£n)n of positive numbers converging to zero, 

^In the case that there are multiple roots, which might be the case for some complex values of 
q, it is still possible to perform a partial fraction decomposition which results in similar but slightly 
different expressions - see Remark 4 in [3]. 
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and, for every n, let {ui)i = and {vj)j = {Vj^^)j be two finite partitions of (0, oo) 

with vanishing mes hi and let (A+)i = {Ap''^)i, {Aj)j = {A- be two finite sets 
of positive weights, shortly to be determined. Set the approximating density equal to 



For a given n, the partitions {ui)i and {vj)j and the weights (Aj^)j, {Aj)j are chosen 
such that the mass of k in the tails and the distance between the target density 
k and the approximating density kn on a closed bounded set not containing zero is 
smaller than that is. 



/ 

Jm.\a 



k{y)dy < En, / {k{y) - kn{y)) dy < En, 

for some open bounded set An,Bn with G Bn C An C M. These requirements can 
be fulfilled since k has finite mass on the set {\y\ > 1} and takes the special form ([1]). 
Set X^") equal to the Levy process with Levy density /c„ and Gaussian coefficient cj^ 

(n) 

= + r] , y\Hy) - kn{y))+dy (9) 

and with drift /i„ chosen such that E[e^''"'^ ^^'^j = e^~'^. The approximating processes 
X^"'^ constructed in this way can be shown to converge weakly to X. H Convergence 
of the barrier option prices corresponding to log price X^") rests on the facl|l| that the 
weak convergence of X*^"-* to X carries over to convergence of the running supremum, 
infimum and the crossing times 

r(")(2;) = inf{t > : X(")(t) < x} and T (x) = mi{t > : X (t) < x} , 

as summarised in the following result. 

Proposition 3.1 Let X be generalised hyper- exponential and let X^"'^ be as defined. 
For T >0, (X(")(r),X^"^(r)) converges in distribution to {X{T),X{T)) as n ^ oo. 
In particular, if X is not of type it holds for t > and z > that 

P(X^'"\t)<z) ^ P(X{t)<z), P{-X^'^Kt) < z) ^ P{-X{t) < z), (W) 
p{T^^){-z) <t) P{T{-z)<t), (11) 



as n ^ oo. 

partition (ti)ilLi of (0, cxd) is an increasing set of times < ti < t2 • • •• The mesh of (ti) is 
defined as maxi<i<m \ti — 

^This convergence takes place in the space _D(R+) endowed with the Skorokhod topology. A proof 
of this fact is given in Appendix [Cl Lemma IC. II 

''a proof is given in Appendix [Cl Lemma IC. 2 1 

''A Levy process is called of type A (see Sato(1999),p. 65) if its Levy measure has finite mass and 
no Brownian motion is present. Type A Levy processes are the compound Poisson processes added 
to a (possibly zero) deterministic drift. 
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4 Laplace transforms of barrier option prices 



In this section we employ the Wiener-Hopf factorisation results to derive the values of 
single barrier options in a hyper-exponential Levy model. We will restrict ourselves 
to down-and-out and down-and-in digital and put options, but related options (such 
as up-and-in call options) can be derived similarly. 

4.1 Digital options 

A digital option is a derivative that pays 1 euro if the price of the underlying asset has 
up or down-crossed a level H before a maturity time T, with pay-off occuring either 
directly at the moment of crossing or at the maturity time T. We will derive the 
values of down-and-out and down-and-in digital options and denote by EDOD and 
EDID the respective prices corresponding to payment at maturity. We also consider 
the American version of the down-and-out digital option (denoted by ADID) that 
pays out at the moment of crossing the barrier. The following result gives the Laplace 
transforms of the digital prices when 5o > H. Denote by / the Laplace transform of 
/: 

Jo 



Proposition 4.1 If q > and H < Sq it holds that 

ADlbiq) = '^^EDlbiq), EDOIDiq) = — EDIDiq). 

q q + r 

where pj = p~{q + r). 

Proof of Proposition 147X1 Since the riskless rate of return is assumed to be constant 
and equal to r, it follows by standard no-arbitrage pricing arguments that the price 
of the down-and-in digital option is given by 

EDID{T) = EDID{T, So) = e-'^P[X{T) < log{H/So)], (12) 

where we switched to log-scale. Also, if we denote h = log{H/So) then 

ADIDiT) = ADIDiT,So) = E[e-''^^''h{xiT)<h}]- 
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Taking now the Laplace transforms of these expressions in T and combining with the 
first-passage distributions given in Section [3] we find if /i < 0: 



.^e 



f'OO f'OO 

EDID{q) = / e-('?+'^)^P[X(T) < h]dT = / e-('?+^)^P[r(/i) < T]dT 
Jo Jo 

1 poo ^ - ^ ^ 

poo r poo 

ADID{q) = / e-'^^E[e'^^(^h{xiT)<H}]dT = E / e-''^-'-^Wl{r(,)<r}dt 
Jo Uo 

= ii?[e-(5+0TW] ^ i+I /■"e-^^{e-^^P[r(/i) < T]}dT. 
Q Q Jo 

The last identity follows by noting that EDOD = e"'"'^ - EDID. □ 
4.2 Knock-out and knock-in put options 

A down-and-out put pays out at maturity T the strike K less the value of the asset St 
(if this difference is positive) with the added feature that the put is worthless if the price 
of the asset has been below a level H by time T. In the next result we derive the Laplace 
transform in the maturity T of down-and-out put DOP{T) = DOP{T; Sq, K, H). The 
value of the down-and-in put is presented in Appendix [Al 

Proposition 4.2 Let q > and write h = log{H/So),£ = log{K/So). Then it holds 
that 

DOPiq) = -^KC^^\l,h) - ^^cW(^,/i), 
q + r q + r 

where, if h < i < 0, 

k~ A;"^ -f- / — N 

and, if h < < i, 



with pf = pf{q + r) and p- = p- [q + r). 
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It is easy to check that, in the case of exponential upward and downward jumps 
(n^ = 1), this result agrees those obtained before by Kou and Wang [26] and Sepp 



Proof of Proposition 14. 2t By the standard theory of no-arbitrage-pricing it follows 
that an arbitrage free price for the down-and-out barrier option is given by 

DOP{T;So,K,H) = e'^^ E[{K - S{T))+1^s(T)>h}\ 
= e-''^KP[X{T) > h, X{T) < I] 

-e-'-^Soi^ie^C^) l{X{T)>h,x(T)<i}], (13) 

where h = \og{H/ Sq) and I = log(K/S'o). Denoting by r = T{q + r) an independent 
random time with an exponential distribution with parameter q + r, the Laplace 
transform DOP{q) of DOP{T; Sq, K, H) can be compactly represented as 

poo 

DOP{q) = K e-^'^^''^^P[X{T) > h,X{T) <£]dT 
Jo 

poo 

-So / e-(''+^)^i^[e^(^)l|x(T)>h,x(T)<n]dr 
Jo 

^ P\X{r) > h, X{t) <l]- ^ii;[e^Ml^^(^)>,_ xir)<i}]m 



q + r q + r 

To calculate the quantities 

C^'\ha) := ^[e^^Ml|^M>h,x(r)<n] 
g{y,b) := E[e''^^^\-^^^^^^^] 

we shall appeal to the following two key properties: 

(a) X{t) and X{t) — XAt) are independent and 

(b) the pairs (X(r),X(r) - X{t)) and (X(r) -X(t),-X(t)) are identically dis- 
tributed. 

In view of (a) it follows that 

e-^^i?[e''(^(-)-^(-))l{^(,)_XM<;.+,}]P(-X(r) e dy) 

— min{A;,0} 

e-'yg{i + y, b)f_xir){y)dy + g{l, b)P[X{T) = 0], 

mm{e,0} 

where we used property (b) in the last line. Inserting the explicit form of the distri- 
bution of Xr, we find that for y > 

9{y,b) = re^^/x(.)(^)dz + P[X(r)=0] 
Jo+ 



£^A^(e(-''^)^-l) + l-£A"- 

i=l "i i=l 
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The next step is to integrate g{£ + y, b) against f-xir)- If ^ < 0, this yields 

k- k+ + 



I _|_ / s 

r^.~i b-pj ■> J-e 

k+ \ k- ..^h 



i=i / j=i 

and in the case that £ > 0, it holds that 



I k^ \ k~ —h 



k+ \ k- f._h 

( n . — K\n ^ 

1=1 I 3=\ " " 

^ + k+ 
Pi 



y j=i J yi=i Pi i=i 

It is then a straightforward matter of some calculus to arrive at the stated expressions. 

□ 

5 Calculating sensitivities 

The analytical formulas for the barrier and digital prices in the previous section suggest 
the possibility of calculating the corresponding sensitivities by direct differentiation. 
In this section we make this idea precise by providing regularity results. The delta 
and gamma sensitivities of the option price with respect to the initial stock price Sq 
will be denoted by 

. dV ^ d'^V 

where V = V{T, Sq) is the value function of the option under consideration (digital or 
down-and-out put). The sensitivity Qy with respect to the maturity T (the theta) is 
defined by a function that satisfies 



V{T) - V{Q) = r ev{s)ds, 
Jo 



for all r > 0. If y is continuously differentiable with respect to T then this definition 

dV 
dT- 



is equivalent to @v = The first result in this direction concerns the theta of a 



digital option: 

Lemma 5.1 Let X be a general Levy process not of type A. Then, for q > 0, Sq > H 

and T>0, 

QEDlDiq) = qEDLDiq), 

QADiDiT) = &EDiDiT) + rEDLD{T). 
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Proof Since the map F : [0, oo) [0, oo) given by F{t) = P{T{x) < t) is an increasing 
function with F{0) = it holds that 

/"OO /"OO 

/ e-''*dF(t) = q / e-'"F{t)dt. (15) 
Jo Jo 

Further, Lemma 49.3 in Sato(1999) imphes that P{T{h) = t) = for /i< 0, t > 0, for 
Levy processes X not of type A, so that F is continuous. If we denote by / a density 
of F it follows that f{q) = qF{q). In view of the form of the pay-off of the digital 
we deduce that &edid{q) = {q + r)EDID{q) — rEDID and the assertion follows. 
The second equation follows from the relation between EDID and ADID given in 
Proposition 14.11 □ 

In the case of a Gaussian component cr > general results are available in the 
literature regarding the smoothness and regularity of the solutions of partial integro- 
differential equations that imply, by the Feynman-Kac representation, that the value 
functions of the barrier options are smooth (see Bensoussan & Lions [7\ and Garroni 
&: Menaldi [21]). More specifically, for a hyper-exponential Levy process with o" > it 
holds that F with F : {t,h) ^ P{X{t) < h) is element of C^'2((0,oo) x (-oo,0)) and 
G : {t,h,k) ^ P{X{i) > h,X{t) < k) satisfies G G (:7i'2'2((o, oo) x (-00,0)^). In view 
of the probabilistic representations ()12|) and ()13p . this result directly implies that the 
delta and gamma of EDID and DOP are well defined. The Laplace transforms of 
the delta and gamma of EDID are given as follows: 

Proposition 5.1 Suppose X is a hyper- exponential Levy process with o" > and let 
q > and H < Sq. It holds that 



^EDID 



TEDiDiq) = -^4i:p-(p--i)^-(5ow^". 

Proof Letting H < c < Sq the fundamental theorem of calculus implies that 



^^-EDID((i) = —J^ J AEDiD{t,y)dydt 

J J AEDiDit,y)dtdy = AEDioiq), 



dSo Jc Jo 

where the change of the order of integration is justified (by Fubini's theorem) since 
Aedid ^ and the last equality follows since EDID{q) = EDID[q; Sq) is contin- 
uously differentiable as a function of 5*0. The form of the delta follows now from 
Proposition 14. 1[ The result for the gamma follows by a similar reasoning as above, 
where the interchange of integration is now justified since EDioit-, So)\ is domi- 

nated by an integrable function (on (0,oo) x {c,Sq)). 

To show the latter fact we note first that since EDID is smooth enough it satisfies 
the partial integro-differential equation that reads in terms of the Greeks of EDID as 

e + ^S'^T + uSA - (A+ + A_ + r)EDID] {t, S) + I{t, S) = 0, 



13 



for t > 0,S > H > with boundary values EDID{t, S) = 1 for S < H,t £ (0, T] and 
EDID{0, S) = 1 if S < H and zero else. Here / is the non-local part of the generator, 
given by 

/oo 
EDID{t, Sey)k{y)dy. 
-oo 

Since EDID < 1 and the measure k{y)dy has finite mass, it thus follows that 

\rEDID{t, S)\ < C(l + \AEDID{t, 5)1 + \eEDID{t, S)\), 

for some constant C. Further, as Aedid > 0, e~'^^\AEDiD{t,y)\dydt is for 

H < c < X equal to AEDinit, y)dy, which is finite in view of the form of the delta. 
Similarly, it follows that e~''*|0(t, y)\dy < oo and the proof is complete. □ 

Given the sensitivities of EDID, those of ADID are calculated as follows: 

Corollary 5.1 Let cr > 0. For q> Q and H < Sq it holds that 

Aadid{T) = AEDiD{T)+r f AEDiD{s)ds, 

Jo 

Tadid{T) = TEDiD{T)+r TEDiD{s)ds. 

Jo 

Proof of Corollary 15. il " The assertions follow the relation between EDID and ADID 
given in Proposition 14.11 and using that f{q)/q is the Laplace transform of f{s)ds. 

□ 

Following an analogous approach, similar formulas can be derived for the sensitivities 
of down-and-out put options; the results are reported in Appendix iBl 

6 Numerical results 

When the stock log-price process is modelled by a generalised hyper-exponential Levy 
process, we use the following algorithm to compute single barrier option prices and 
corresponding sensitivities: 

Algorithm 

1. Approximate the target Levy density by a hyper-exponential density. 

2. Calculate the Laplace transforms of prices and sensitivities using 
the formulas in Sections |4] and [5]. 

3. Apply a Laplace inversion algorithm to the result of step 2. 
Comments 

Ad 1. Similarly as in Asmussen et al. [5] we specified the form of the density ([5]) by 
fixing the number of exponentials and the mean jump sizes {af)~^, {aj)~^ and then 
using a least squares algorithm to determine the values of the remaining parameters 
pf, pj and A^, \J that minimize the squared distance between the density (jS]) and 



14 



the target density. To improve the accuracy of the fit one could employ methods from 
the approximation theory of real valued functions, or alternatively follow the approach 
developed by Feldmann and Whitt [19]. 

Ad 2. The roots of the characteristic exponent that appear in the formulas of 
the Laplace transforms were calculated using Laguerre's method (see e.g. Numerical 
Recipes [32]). 

Ad 3. The quantity of interest V can be expressed in terms of the Laplace transform 
V{q) by the Bromwich integral 



-1 /-c+ioo 

ViT) = — / e'>'^V{q)dq 



(16) 



To evaluate the integral (I16p we employed the algorithm investigated by Abate and 
Whitt [T], which is based on approximation of the integral (|16p by Euler summation 
with the trapezoidal rule that leads to the series 

M 



where 



Sn{T) 



-Re 



oA/2 



k=l 



n^) 



Abate and Whitt [Ij recommend the values M = 15, = 11, A = 18.4, where 
should be increased if necessary. See [U p. 38-39] for a discussion about the role these 
parameters play in controlling the error. 

We have implemented the above algorithm in the context of a variance gamma 
(VG) model and a normal inverse gaussian (NIG) model to calculate the prices and the 
sensitivities of digital barrier options and down-and-out put options. The parameters 
of the VG and NIG models were determined by calibrating the models to market 
quoted Stoxx50E call prices on 16 June 2006, using the FFT algorithm proposed by 
Carr and Madan [TS]. This procedure yielded the following parameter values C = 
0.925, G = 4.667 and M = 11.876 for the variance gamma model and a = 8.858, 
/3 = —5.808, 6 = 0.174 for the normal inverse Gaussian model. We employed these 
parameters in the subsequent calculations of barrier option prices and sensitivities. 

As an initial step, the Levy densities of the VG and NIG processes were approxi- 
mated by hyper-exponential densities of the form ^ with 14 different exponentials, 7 
upward phases (n"^ = 7) and 7 downward phases (n~ = 7). Note that the parameters 
af and aj are equal to the reciprocals of the mean sizes of upward and downward 
jumps of the (hyper-exponential) log-price, respectively. The values of af and aJ were 
fixed as in Tabled] The remaining parameters A"*", A~, and the pf, pj were determined 
by using a least squares algorithm to minimize the sums of squares of the difference 
between the target density and the hyper-exponential density ([5|) over a time-grid 
inside the interval [0.1,5]. The resulting parameter values of the hyper-exponential 
density are specified in Table [T] 

We simulated paths of variance gamma and normal inverse Gaussian processes by 
employing the representations of these processes as random time changes of a Brownian 
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Table 1: Parameter values for the approximation of the VG (C = 0.925, G 
4.667, M = 11.876) and NIG (a = 8.858,/? = -5.808, <5 = 0.174) Levy densities 
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Figure 1: Prices and sensitivities of a down-and-out put option with strike K — 3500 EUR 
and barrier H set at 60% of K in an exponential NIG model. Price (left up) - Delta (right 
up), Gamma (left-down) and Theta (right down). Results from the transform algorithm are 
indicated with a x and Monte Carlo results with a o. 



motion by an independent gamma and an independent inverse Gaussian subordinator, 
respectively. The sensitivities were calculated by approximating the derivative by a 
finite difference and subsequently evaluating this using Monte Carlo simulation. 

Tables [5] and E] and Figures [T] and [5] report prices and greeks for down-and-in 
digital and down-and-out put options calculated by using Algorithm 1 and Monte- 
Carlo simulation. The options are priced at different spot levels that are expressed 
as a percentage of the spot price on June 16th 2006 (S'o = 3500). The value of the 
sensitivities is expressed as a fraction of S'o- Down-and-in digital options are priced 
under the VG model with barrier level H set at 60%, and down-and-out puts under 
the NIG model with a barrier H and a strike K set respectively at 60% and 100%. 

The results show a general agreement between the Monte Carlo simulation results 
and those from the semi-analytical approximation. The relative errors for the prices 
are less than 0.6% for VG and about 2% for NIG, and for the delta about 2% - 2.5% for 
both, if, for NIG, we stay some distance away from the barrier. We also observe that. 
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(157, 
(141, 



166) 
149) 
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-5.32 
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-2.65 
-1.93 
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0.10 
0.33 
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0.43 
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0.35 
0.30 
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0.12 
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(-9.15, 

(-5.47, 
(-3.75, 
(-2.65, 
(-2.06, 
(-1.35, 
(-0.95, 
(-0.54, 
(-0.19, 
(0.04, 
(0.30, 
(0.37, 
(0.57, 
(0.55, 
(0.51, 
(0.52, 
(0.44, 
(0.43, 
(0.36, 
(0.29, 
(0.25, 
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(0.14, 
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-8.98) 
-5.31) 
-3.60) 
-2.51) 
-1.92) 
-1.22) 
-0.91) 
-0.50) 
-0.14) 
0.12) 
0.38) 
0.49) 
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0.67) 
0.63) 
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0.12) 
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-360 
-373 
-369 
-354 
-333 
-309 
-281 
-250 
-217 
-181 
-145 
-110 
-77.5 
-49.2 
-25.6 
-6.7 
8.0 
18.9 
26.9 
32.4 
36.2 
38.4 
39.7 
40.1 
39.9 
39.3 
38.4 
37.4 
36.1 
34.8 



(-365,-342) 
(-377,-356) 
(-373,-351) 
(-358,-338) 
(-338,-318) 
(-316,-297) 
(-286,-267) 
(-255,-237) 
(-224,-206) 
(-187,-169) 
(-147,-129) 
(-112,-95) 
(-80,-63) 
(-54,-38) 
(-30,-14) 
(-12,3) 
(3,18) 
(15,30) 
(21,35) 
(25,39) 
(30,43) 
(33,45) 
(34,46) 
(34,46) 
(32,44) 
(31,42) 
(30,40) 
(29,38) 
(28,37) 
(27,36) 



Table 2: The prices and sensitivities of a down-and-out put in an exponential NIG Levy model 
with maturity T = 1 year, strike K = 3500 and barrier H set at 60% of K. The interest and 
dividend rates are taken to be r = 0.03 and d = 0. All columns except that of the price are 
expressed in percentage figures of K. The columns with LT contain the results obtained using 
the transform algorithm, whereas MC refers to Monte Carlo results, which are reported in the 
form of a 95% confidence interval. The columns with RE contain the relative errors, which 
is computed as |MC^ — LT\/LT respectively, where MC^ is the mid-point of the confidence 
interval. 
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Figure 2: Prices and sensitivities of an American down-and-in digital option with a barrier H 
set at 60% of EUR 3500 in an exponential VG model. Price (left up) - Delta (right up), Gamma 
(left-down) and Theta (right down). Results from the transform algorithm are indicated with 
a X and Monte Carlo results with a o. 



close to the barrier, the result for the gamma falls outside the confidence interval. In 
the case of the NIG Levy model, the larger errors close to the barrier may be explained 
by the fact that the option price is not smooth under this model whereas it is for the 
approximating hyper-exponential model (see [13] for an analysis of the smoothness of 
option prices under an NIG or, more generally, a regular Levy process). 

The convergence of Monte Carlo estimators of quantities involving first passage 
is known to be very slow, requiring a small mesh of the time grid. Since this effect 
is magnified when calculating the numerical derivatives, a large number of paths and 
time-steps was needed to obtain convergence. For the Monte Carlo calculations we 
used 1,000,000 paths and 20,000 time steps per year, which made the calculation of 
the greeks using Monte Carlo computationally demanding and time-consuming. In 
a C-|--|- programme on a 3189 Mhz computer the valuation of all option prices and 
sensitivities took 40 seconds for a digital option and 55 seconds for a down-and-out 
put option, compared to a couple of hours when using Monte Carlo. Here we should 
remark that we did not employ any special methods that could considerably have 
improved the speed of convergence of the Monte Carlo simulation, such as Malliavin 
weights, likelihood ratio methods or sampling techniques such as the bridge (see e.g. 
Glasserman [23]). We noted that if the spot price was close to the barrier the Monte 
Carlo estimators of the greeks tended to be unstable. 



7 Convergence of prices and sensitivities 

A generalised hyper-exponential Levy process can, as we have shown in Section 13.21 
be approximated arbitrarily closely by a hyper-exponential jump-diffusion, by speci- 
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(-14.4, 
(-12.5, 
(-10.8, 
-9.4, 
-8.4, 
-7.4, 
-6.5, 
-5.8, 
-5.1, 
-4.6, 
-4.1, 
-3.7, 
-3.3, 
(-3.0, 
-2.7, 
-2.4, 
-2.2, 
-2.0, 
-1.7, 
-1.6, 
-1.5, 



-106) 
-81.0 
-64.3 
-52.1 
-43.0 
-35.8 
-30.1 
-25.4 
-21.6 
-18.4 
-16.0 
-13.8 
-12.0 
-10.4 
-9.1 
-8.0' 
-7.1 
-6.2 
-5.5 
-4.9' 
-4.4' 
-3.9^ 
-3.5 
-3.1 
-2.8 
-2.5' 
-2.2 
-2.1 
-1.9' 
-1.7' 
-1 
-1.4 



2.9 
1.9 
1.3 
0.7 
0.3 
0.0 
0.4 
0.4 
0.8 
0.9 
0.2 
0.1 
0.5 
0.4 
0.6 
0.1 
0.7 
0.6 
0.6 
1.0 
0.3 
0.1 
1.4 
1.4 
0.3 
0.2 
1.8 
0.2 
0.5 
0.6 
2.3 
2.2 



43.2 
27.4 
19.4 
14.6 
11.4 
9.06 
7.34 
6.02 
4.97 
4.13 
3.46 
2.90 
2.45 
2.08 
1.77 
1.52 
1.30 
1.12 
0.97 
0.84 
0.73 
0.64 
0.56 
0.49 
0.43 
0.38 
0.34 
0.30 
0.27 
0.24 
0.21 
0.19 



(44.4 
(27.4 
(19.5 
(14.1 
(10.9 
(8.55 
(7.07 
(5.35 
(4.85 
(3.22 
(3.05 
(2.36 
(2.14 
(2.06 
(1.26 
(1.28 
(1.08 
(0.92 
(0.64 
(0.54 
(0.47. 
(0.40 
(0.31 
(0.28 
(0.25 
(0.23 
(0.22 
(0.21 
(0.17. 
(0.16 
(0.15 
(0.14, 



45.4) 
28.3) 
20.9) 
15.2) 
12.2) 
9.57) 
7.99) 
6.17) 
5.89) 
4.15) 
3.89) 
3.11) 
2.82) 
2.64) 
1.81) 
1.79) 
1.38) 
1.29) 
0.99) 
0.89) 
0.84) 
0.78) 
0.72) 
0.63) 
0.61) 
0.48) 
0.41) 
0.37) 
0.33) 
0.30) 
0.26) 
0.23) 



19.0 
18.5 
17.6 

16.6 
15.4 
14.2 
13.1 
12.1 
11.1 
10.2 
9.34 
8.57 
7.87 
7.22 
6.64 
6.11 
5.63 
5.19 
4.79 
4.42 
4.09 
3.78 
3.51 
3.25 
3.02 
2.81 
2.62 
2.44 
2.27 
2.12 
1.98 
1.86 



(16.3,19.4) 
(16.3,19.3) 
(16.2,19.1) 
(14.8,17.6) 
(13.2,15.8) 
(12.5,15.0) 
(11.6,14.0) 
(11.2,13.5) 
(10.1,12.2) 
(9.4,11.4) 
(8.2,10.2) 
(7.6,9.4) 
(6.7,8.4) 
(6.6,8.2) 
(6.2,7.8) 
(5.9,7.4) 
(5.4,6.8) 
(4.9,6.3) 
(4.3,5.2) 
(4.1,5.3) 
(3.6,4.8) 
(3.5,4.6) 
(3.2,4.3) 
(3.0,4.0) 
(2.7,3.7) 
(2.6,3.4) 
(2.3,3.2) 
(2.2,3.1) 
(2.0,2.9) 
(1.8,2.6) 
(1.7,2.6) 
(1.6,2.4) 



Table 3: The prices and sensitivities of an American down-and-in digital option with maturity 
T = 1 year barrier and H set at 60% of 3500 in an exponential VG Levy model. The interest 
and dividend rates are taken to be r = 0.03 and d = 0. All columns except that of the prices 
are expressed in percentage figures of 3500. The columns with LT contain the results obtained 
using the transform algorithm, whereas MC refers to Monte Carlo results, which are reported 
in the form of a 95% confidence interval. The columns RE contain the relative error, which 
is computed as \MC^ — LT\/LT respectively, where MC^ is the mid-point of the confidence 
interval. 
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fying appropriate values for the parameters of the latter process. Under such a model 
analytical expressions were derived in Sections S] and [5] for the values EDID, EDOD, 
ADID, DOP of digital and the down-and-out put options, and the corresponding sen- 
sitivities. The results in Section [6] provided numerical evidence to show that a good 
approximation of prices and sensitivities in a GHE Levy model can be obtained by 
carrying out the computations in a suitably chosen approximating hyper-exponential 
jump-diffusion model. In this section we provide further justification for the algo- 
rithm by proving a convergence result, that is, for a given GHE Levy process X and 
a sequence of HEJD processes weakly converging to X (as constructed in Sec- 
tion [32]), we will show that the corresponding prices EDID'~'^\ ED0D'-'^\ ADID^'^\ 
DOP^"'^ of digital and down-and-out put options converge pointwise to those in the 
limiting model X, as n — > oo: 

Proposition 7.1 Suppose that X is GHE and is not of type A and let V be any of 
EDID, ADID, EDOD, DOP with V^^^ the corresponding approximation. Then for 
Sq > H,T > 0, it holds that 

y(")(r, So) V{T, So) as n ^ oo. (17) 

Proof of Proposition \ 7.1\ The convergence of EDID^-^\T,So) to EDID{T,So) IS a 



direct consequence of the convergence of p(r(")(/i) < t) (Proposition El]). Further, 
by interchanging the order of integration it follows that the value function ADID can 
be written as 

r-oo 

ADID{T) = / re^''*P{T{h) < T A t)dt. 
Jo 

Thus, in view of the dominated convergence theorem and again Proposition 13.11 it 
follows that ADID^'^\T,So) ADID{T,So). Also, since {K - e^)+ is a contin- 
uous bounded function. Lemma [C.21 implies that E[{K — e^'^"^^'^'^)^'i-{T'."'){h)<T}] ~ 
DOP^^-^T; So) converges to DOP{T; So) as n ^ oo. □ 

Note that Proposition 17.11 implies that the finite difference approximations of the 
sensitivities of the digital and down-and-out put options also converge pointwise as 
n — > oo for any spot price away from the barrier H and any maturity T > 0. To 
rigorously prove pointwise convergence of the sensitivities, uniform estimates would 
be needed of the errors in p7|) . It is worth noting that, while for Levy processes 
with positive Gaussian component first-passage probabilities and value-functions of 
barrier options are smooth up to the barrier, such is not necessarily the case for a 
Levy process without Gaussian component, especially at the barrier where the spatial 
derivatives may be infinite. However, when suitably smoothed, the sensitivities do 
converge pointwise: 

Corollary 7.1 For some e > let V''{s) = f V(T,y)(p^{s — y)dy, where (f)^ is a 
probability density with support in (— e, e). Then it holds that 

^y(r^){So,T) Av,{So,T), T^(„){So,T) TvSSo,T), 
for So > H + e,T > 0, as n ^ oo. 
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A similar statement holds true for the theta, replacing smoothing in space by smooth- 
ing in time. Note that, if s ViT, s) is on {H, oo), then Av;(5o, T) — Ay(S'o, T), 
Sq > H, can be made arbitrarily small by choosing e small enough (similar remarks 
apply to the other greeks) 

Proof In view of the definition of v}^^ and integration by parts it follows that 

A^(„)(5o,T) = AyW(5„,r) = I V^^\ymSo - y)dy. (18) 

By the dominated convergence theorem and Proposition 17.11 it thus follows that the 
rhs of ([TH]) converges to J V{y)(j)'^{So — y)dy = Ay^So, T) as n ^ oo. The convergence 
of the gamma follows by a similar reasoning. □ 

8 Conclusion 

In this paper we developed an efficient algorithm to compute prices and sensitivities of 
barrier options driven by an exponential Levy model of generalised hyper-exponential 
type. We showed that the latter class contains many of the Levy models that are 
employed in mathematical finance. We first approximated the target Levy measure 
by a hyper-exponential one. Subsequently, for log-price processes in this class, jump- 
diffusions with hyper-exponential jumps, we derived analytical expressions for the 
prices and sensitivities (greeks) of digital, knock- in and knock-out option prices, up 
to a single Laplace transform. Inversion of this Laplace transform yielded fast and 
accurate results for the option prices and sensitivities. We proved convergence of this 
algorithm. To provide a numerical illustration we implemented the algorithm for the 
VG and NIG Levy models, approximating the NIG and VG Levy densities by hyper- 
exponential ones with 7 upward and 7 downward phases. Compared with Monte Carlo 
simulation results we found relative errors of about 0.5-2.5% for prices and deltas some 
distance away from the barrier. What the rate of convergence of this algorithm is when 
increasing the number of terms in the approximation, and how this rate depends on 
the quality of the approximation of the target density, the parameters and the distance 
to the barrier, are open questions that are left for future research. 
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APPENDIX 



A Additional pricing formulas 
A.l Upward digital options 

The arbitrage free prices of the corresponding up-and-in digitals (contracts that pays 
out £1 if an upper barrier is crossed, either at T or at the moment T^{h) of up- 
crossing) are given by 

EUID{T,So,H) = e-'^^P[5(T) > H], AUID{T,So,H) = S[e-'^^''Wl^5(^)>^^], 

where ^(r) = supo<^<y /S'(s) denotes the running supremum of S up to T. Their 
respective Laplace transforms in T follow by applying the formulas of the down-and- 
in digital options to the process —X: 



1 

EUlb{q) = Afe-P^^, AUID{q) = ^^^EUlb{q), h>0. 



i=l 



A. 2 Down-and-in put 

The Laplace transform DIP{q) of the arbitrage free price of the down-and-in put 
option DIP{T; Sq, K, H) with strike K and barrier H can be decomposed as 



q + r q + r 

where h = log{H/So), k = log{K/So) and D^^\h, k) is defined as 

D(-''\h,k) := E[e'"'^^h{x(r)<h,X(r)<k}]- 

Reasoning as in Proposition 14.21 we find that 

D^'\h,k) = / e-^^ii;[e^(^W-^«)l|^(,)_XM<fc+,}]^'(-^(r)Gdy). 



-h 

After some calculations we arrive at 

fc- k+ 



f^,t^, {Pj - Pt)ib - Pi) 



22 



B Sensitivities for down-and-out-puts 

Proposition B.l Suppose X is a hyper- exponential Levy process with a > and let 
q > and H < min{i^, Sq}. 

(i) It holds that QDOp{q) = qDOP{q) - {K - Sq)+ . 

(a) If So < K, it holds that 



where 



k- 



k+ 



1=1 



(p; -pt)ib-pt) 



+ 



Pj){b-Pt) (Pj -pT)ib-pJ) 



+ 



1 



b-p, 

(Hi) If K < So, it holds that 



^DOp{q) 



Y.^{p]fS^~"\KFf-Ff^] 



i=l 



+ Y.^tiptfipt - l)S'J~\KGf - Gf], 



1=1 



where 64 



1 - A^T^ and 5. = l- E)=i AJ 



and 



b- Pi ^ 



i=l 



K^-Pt 
b^ 



Pt 



j=i ' Pj - pt 



{-pl)K^-pt Rf^ 
ipj - Pt)ib - pf) 

1 



H^-f^ b 



Pj)ib-Pt) 



Proof of Proposition \B.1[ The formula of the Laplace transform of the theta is a 
consequence of integration by parts and the fact that h, k) is continuously differ- 
entiable on (0, 00). Those of the Laplace transforms of the delta and the gamma follow 
by differentiating the expressions in Proposition 14.21 following a reasoning analogous 
to the one in Proposition 15.11 □ 
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C Proofs of the weak convergence 



Let X be a generalised hyper-exponential Levy process and define X^") as in Section[3l 
Refer to Jacod and Shiryaev [2^ Ch. VI] for background on the Skorokhod topology 
and weak convergence of stochastic processes. Let Z)(M+) denote the space of real 
valued right-continuous functions with left-limits on 

Lemma C.l X^"^ converge weakly to X in the Skorokhod topology on D{M^), as 
n — > oo. 

Proof . With regard to Corollary VIL3.6 in Jacod and Shiryaev [241 Ch. VI] it follows 
that it is equivalent to show that converges in distribution to X{1), and that 

the latter assertion follows if it holds that 



Hn^lJ', J g{u)kn{u)du ^ J g{u)k{u)du, 
crn~^ J x'^kn{x)dx ~^ ^"^ ^ J x'^k{x)dx, 

for any continuous bounded function g that is zero around zero, as n — > oo. As a 
consequence of the definitions of kn and cr^ given in ([8]) and ([9]) it is not hard to verify 
that these three conditions are indeed satisfied. □ 

Lemma C.2 (X(")(r), X^"\r)) converges in distribution to {X{T),X{T)). 

Proof . It follows from Proposition VI. 2. 4 and remark VI. 2. 3 in Jacod and Shiryaev 
[24j that for each fixed T such that Auj{T) = 0, the function / : -D(M+) M given by 

f : uj ^ (cj(T), sup{max{u;(s), 0}, s < T) 

is continuous (in the Skorokhod topology). In particular, if w^"^ ^ u; in the Skorokhod 
topology, then f{uj^'^'>) /(w), as n ^ oo. In view of the fact that Levy processes are 
continuous at each fixed time a.s., it holds that P(AX(T) 7^ 0) = 0, and the asserted 
convergence follows in view of Lemma IC.ll □ 
Proof of Proposition \3.1\ The convergence of (X(")(r),X^"^(T)) is proved in Lemma 
IC.2I If X is not of type A, it is shown in Lemma 49.3 in Sato(1999) that X{T) and 
—X{T) are continuous random variables on (0, 00) (with possibly an atom in zero if 
X is of bounded variation). The statement (jlOp thus follows in view of the continuity 
theorem. Since P{T{x) < T) = P{X{T) < x), also 1^ follows. □ 
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